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ABSTRACT 

A detailed analysis of a coronal loop oscillation event is presented, using data 
from the Atmospheric Imaging Assembly (AIA) onboard the Solar Dynamics Ob- 
servatory (SDO) for the first time. The loop oscillation event occurred on 2010 
Oct 16, 19:05-19:35 UT, was triggered by an M2.9 GOES-class flare, located in- 
side a highly inclined cone of a narrow-angle CME. This oscillation event had a 
number of unusual features: (i) Excitation of kink-mode oscillations in vertical 
polarization (in the loop plane); (ii) Coupled cross- sectional and density oscil- 
lations with identical periods; (iii) no detectable kink amplitude damping over 
the observed duration of four kink-mode periods [P = 6.3 min); (iv) multi-loop 
oscillations with slightly (~ 10%) different periods; and (v) a relatively cool loop 
temperature of T ^ 0.5 MK. We employ a novel method of deriving the elec- 
tron density ratio external and internal to the oscillating loop from the ratio 
of Alfvenic speeds deduced from the flare trigger delay and the kink-mode pe- 
riod, i.e., Ue/ui = (va/vac)^ = 0.08 ± 0.01. The coupling of the kink mode and 
cross-sectional oscillations can be explained as a consequence of the loop length 
variation in the vertical polarization mode. We determine the exact footpoint 
locations and loop length with stereoscopic triangulation using STEREO/EUVI- 
A data. We model the magnetic field in the oscillating loop using HMI/SDO 
magnetogram data and a potential field model and find agreement with the seis- 
mological value of the magnetic field, -Bfcmfc = 4.0 ± 0.7 G, within a factor of 
two. 



Subject headings: Sun: Flares — Sun : Corona — Sun: Extreme Ultra- Violet 
(EUV) — Sun : oscillations — Sun : waves 
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1. INTRODUCTION 

Propagating waves and standing waves (eigen-modes) in coronal plasma structures be- 
came an important tool to probe the physical parameters, the dynamics, and the magnetic 
field in the corona, in flare sites, and in coronal mass ejections (CMEs). Recent reviews 
on the theory and observations of coronal seismology can be found in Roberts and Nakari- 
akov (2003), Erdelyi et al. (2003), Roberts (2004), Aschwanden (2004, 2006), Wang (2004), 
Nakariakov and Verwichte (2005), Banerjee et al. (2007), Andries et al. (2009), Ruderman 
and Erdelyi (2009), and Taroyan and Erdelyi (2009). Substantial progress was accomplished 
in applying MHD wave theory to coronal observations with previous instruments, such as 
the discovery of global waves with EIT/SOHO (Thompson et al. 1998, 1999), fast kink-mode 
loop oscillations with TRACE (Aschwanden et al. 1999; Nakariakov et al. 1999), fast sausage 
mode oscillations in radio wavelengths (Roberts et al. 1984; Asai et al. 2001; Melnikov et 
al. 2002; Aschwanden et al. 2004), slow (acoustic) mode oscillations with SUMER/SOHO 
(Wang et al. 2002; Kliem et al. 2002), slow (acoustic) propagating waves with UVCS/SOHO 
(Ofman et al. 1997) and Hinode (Erdelyi and Taroyan 2008), EIT/SOHO (DeForest and 
Gurman 1998), and TRACE (De Moortel et al. 2002a,b), fast Alfvenic waves with SECIS 
(Williams et al. 2001; Katsiyannis et al. 2003), or fast kink waves with TRACE (Verwichte et 
al. 2004; Tomczyk et al. 2007). Sausage oscillations observed and identified directly in cross- 
sectional area change of a solar magnetic fiux tube are reported by Morton et al. (2001). 
However, temporal cadence of space-borne EUV imagers (such as EIT/SOHO, TRACE, 
STEREO) was mostly in the order of 1-2 minutes, which is just at the limit to resolve fast 
MHD mode oscillations (with typical periods of 3-5 minutes) and is definitely too slow to 
resolve or even detect fast MHD waves that propagate with Alfvenic speed. An Alfven wave 
with a typical coronal speed of va ~ 1000 km s~^ traverses an active region in about 1 minute. 
With the advent of the Atmospheric Imaging Assembly (AIA) onboard the Solar Dynamics 
Observatory (SDO), which provides a permanent cadence of 12 s, we have an unprecedented 
opportunity to study the exact timing of the excitation mechanisms of coronal MHD waves 
and oscillations, which often involve an initial impulsive pressure perturbation in a fiare and 
CME source site, that launches various MHD waves and oscillations in surrounding resonant 
coronal structures (loops, fans, CME cones, and cavities). 

Here we conduct a first AIA study on a coronal loop oscillation event, observed on 2010 
October 16, which exhibits a favorable geometry, unobstructed view, prominent undamped 
oscillations, an unusual coupling of kink-mode and cross-sectional (and density) oscillations 
(not noticed earlier), and a rare case of vertical kink-mode polarization. In Section 2 we 
present various aspects of the data analysis and modeling, while theoretical and interpreta- 
tional aspects are discussed in Section 3, with the major findings and conclusions summarized 
in Section 4. 
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2. DATA ANALYSIS 
2.1. Instrument 

The Atmospheric Imaging Assembly (AIA ) instrument onboard the Solar Dynamics Ob- 
servatory (SDO) started observations on 2010 March 29 and produced since then continuous 
data of the full Sun with four 4096 x 4096 detectors with a pixel size of 0.6", corresponding 
to an effective spatial resolution of ~ 1.6". AIA contains 10 different wavelength channels, 
three in white light and UV, and seven EUV channels, whereof six wavelengths (131, 171, 
193, 211, 335, 94 A) are centered on strong iron lines (Fe VIll, IX, Xll, XIV, XVI, XVlll), 
covering the coronal range from T ^ 0.6 MK to > 16 MK. AIA records a full set of near- 
simultaneous images in each temperature filter with a fixed cadence of 12 s. Instrumental 
descriptions can be found in Lemen et al. (2011) and Boerner et al. (2011). 

2.2. Observations and Location 

A major fiare of GOES class M2.9 occurred on 2010 Oct 16, 19:07-19:12 UT at location 
W26/S20 ("1-390" west and —410" south of Sun center), which triggered a number of loop 
oscillations in the westward direction of the active region (NOAA 1112). In this study we 
focus on the detailed analysis of a loop at the loop apex position -1-698" west and —243" south 
of Sun center, which displays prominent oscillations. The location of this loop with respect 
to the flare center, is shown in Fig. 1. The oscillating loop is discernible as a faint semi- 
circular structure in the logarithmically-scaled intensity image in 171 A (Fig. 1, top panel), 
or even clearer in the difference image (19:22:36 UT - 19:21:00 UT) in Fig. 1 (bottom panel), 
where the times were chosen at the maximum and subsequent minimum of an oscillation 
period. The 171 A intensity image (Fig. 1) shows also irregular moss-like structure in the 
background of the oscillating loop of similar brightness, which poses some challenge for exact 
measurements of the loop oscillation parameters, because the background is time-variable, 
even on the time scale of the oscillation period. 

2.3. Transverse Loop Oscillations 

Loop oscillations are traditionally investigated most easily in time-difference movies. 
[Movies of this oscillation event in 171 A intensity and running- difference format are available 
as supplementary data in the electronic version of this journal]. However, a variety of time- 
differencing schemes can be applied in order to enhance the best contrast. We explore a 
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variety of time-differencing schemes in Fig. 2, for a data stripe oriented perpendicular to the 
loop axis at its apex with a length of 30 pixels and a width of 10 pixels (indicated with a 
small rectangle in Fig. 1). We construct time-slice plots with Ut = 150 time frames on the 
X-axis (covering the time interval from 19:05 UT to 19:35 UT with a cadence of At = 12 
s) and a spatial dimension in direction transverse to the loop on the y-axis (with Uy = 30 
pixels), averaged over the = 10 pixels of the stripe width (parallel to the loop). We show 
5 different differencing schemes of this time shce in Fig. 2, using a high-pass filter (Fig. 2 top 
panel), a baseline difference (Fig. 2, second panel), and a one-sided (Fig. 2, third panel), a 
symmetric (Fig. 2, fourth panel), and a running- minimum difference scheme (Fig. 2, bottom 
panel), which is defined as 



so it subtracts a running minimum evaluated within a time interval with a length of 2k pixels 
symmetrically placed around every time slice. Each method has its merits and disadvantages, 
as it can be seen in Fig. 2. The biggest challenge is the non-uniformity and time variability of 
the background. An additional complication is the presence of fainter secondary oscillating 
loops, which appear hke "echoes" in the time-shce plots. For further analysis we adopt 
the running-minimum differencing scheme (Fig. 2, bottom), which appears to have the best 
signal-to-noise ratio of the oscillating features. 

The measurement of the loop oscillation amplitude variation a{t) as a function of time t 
can be done (i) by localizing the cross-sectional fiux maxima in running- difference time-slice 
plots, (ii) by cross-correlation of subsequent time slices, or (iii) by fitting a Gaussian profile 
to the cross-sectional fiux profiles. We find that the first and the latter method are most 
robust. From the running- minimum time-slice plot (Fig. 3 top frame) we perform fits of 
Gaussian profiles Ffit{s,t) to the observed cross-sectional fiux profiles F{s,t) in each time 
slice t (using the standard GAUSSFIT.PRO routine in the IDL software). 



which yields the four coefficients of the peak fiux f{t), the oscillation amplitude a{t), Gaus- 
sian width (7s (t), and mean background flux b{t) for each time t. The 4-parameter flts The 
cross-sectional flux proflles F{s,t) and the Gaussian flts Ffit{s,t) are shown in Fig. 4 for 
each time in the interval between ti=19:05 UT and ti5o=19:35 UT, while a corresponding 
time-shce with Gaussian flts is rendered in color scale in Fig. 3 (second panel). The average 
Gaussian loop width during the oscfllation period is a g — 2.1 Mm, which corresponds to a 
FWHM loop width of w = (7,2V21n2 = 4.9 ± 0.6 Mm. 

We are fltting now a sinusoidal function with a linear drift to the location of the oscilla- 
tion amplitudes a{t) (crosses in Fig. 3, fourth panel), using the Powell optimization routine 



AF(ti, yj) = F{ti, yj) - min[F{ti_k, Vj), F{ti+k, Vj)], 



(1) 
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(Press et al. 1986) from the IDL software, 



afit{t) = ao + oisin 



( 



2n{t - to) 
P 



) 



+ 02 



{t - to) 
P 



(3) 



for which we find a midpoint position oq = 6.3 Mm, a drift velocity a2/P = 0.8 km/s, an 
oscillation period of P = 395 s (6.4 min), an oscillation amplitude oi = 1.8 Mm, and a 
sinusoidal onset time of to — 393 s after the start of the time slice at 19:05:00 UT, i.e., at 
19:11:33 UT. The onset time of the oscillation will be important to measure the exciter speed 
of the trigger. The fit of the sinusoidal amplitude function a/it(i) to the measured amplitude 
a{t) is shown in Fig. 3 (fourth panel). The fitted function with a constant amplitude ai 
appears to be appropriate for the duration of Npuise = (^iso ~ t?,'i)/P = 1407/395 = 3.6 
oscillation periods, since we do not observe any significant damping of the amplitude during 
this time interval. 



The projected loop shape is close to a semi-circular geometry (Fig. 1, bottom), and thus 
we can assume that the loop plane is near the plane-of-sky or nearly perpendicular to the 
line-of-sight. The location of the loop curvature center is at a distance of ~ 740" from Sun 
center or 0.77 solar radii, which corresponds to a heliographic angle of a = 50° from disk 
center. 

The full 3D geometry of the loop can be obtained from the combination of the EUVI 
instrument onboard STEREO and AIA observations, a procedure that we carry out for the 
first time here. The loop was in the field-of-view of STEREO/A(head) at this time, but 
was occulted for STEREO/B. The STEREO/A spacecraft was located on 2010 Oct 16 at 
a separation angle of aA — 83.583° to the east of Earth, at a latitude of — —0.119° 
from the Earth echptic plane. In Fig. 5 we show nearly contemporaneous AIA and EUVI/ A 
difference images of the loop, which clearly show the oscillatory motion of the loop, after 
highpass- filtering of the EUVI/ A image. Unfortunately, EUVI/ A observed only in a different 
wavelength of 195 A at this time, while the oscillation is best visible in the 171 A channel in 
AIA. EUVI/A had also a lower cadence (?^ 5 min vs. 12 s in AIA) and the spatial resolution 
of EUVI (1.6" pixels) is about three times coarser than AIA (0.6" pixels). Nevertheless, the 
image quality is sufficient to approximately determine the 3D loop geometry. We subtracted 
the earlier image from the later image, and thus a density increase in the difference images 
(white in Fig. 5) indicates an inward loop motion (in the AIA image) and a correlated density 
compression (in the EUVI/A image). We rotate the 2D-coordinates of the loop traced in AIA 
(Fig. 5 left) into the coordinate system of EUVI/A with variable heights and inchnation angle 
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of the loop plane. By matching the position and direction of the loop ridge in EUVI/A we 
obtain the absolute height range of the traced loop segment, i.e., 21.7 < hsegm < 37.4 Mm. In 
order to locate the positions of the footpoints we extrapolate the traced loop segment in both 
directions and define the positions of the loop footpoints where the coplanar extrapolation 
intersects with a height h = above the solar surface. The so-defined extrapolated footpoint 
positions are found at Fi = (685", —305") (south of traced loop) and F2 = (615", —268") 
(east of trace loop) with respect to Sun center (Fig. 5). The apex or midpoint of the traced 
loop segment (at s — Lioop/2) is located at position {x apex ■,y apex) — (698", —243"), for which 
we show time-slice plots of the oscillation in Figs. 2-4 (i.e., segment #6 in Fig. 6). The apex 
location will also be used to define the arrival time of the exciting wave and starting time 
of the loop oscillation in Section 2.11. The inclination angle of the loop plane to the local 
vertical is found to be ^9 ^ 20° ± 20°, but cannot be determined more accurately because of 
the short loop segment detectable in EUVI/A. 

From the absolute 3D coordinates (xj, 2^1), i = 1, n of the stereoscopically triangu- 
lated loop we can calculate the full loop length L^oop, 

n— 1 

Lioop = ^ - + {yi+i - + {zi+i - z^Y] , (4) 

i=0 

for which we find L^op = 163 Mm. The traced loop segment over which amplitude oscillations 
are clearly visible covers the fraction from Si/ L = 0.23 to S2/L = 0.78 of the total loop length 
and has only a length of Lgegm = 123 Mm. If we approximate the 3D loop geometry with a 
coplanar semi-circular shape, the loop curvature radius is estimated to be rioop ~ Lioop/'n = 
52 Mm. 

The plane of transverse loop oscillations with respect to the average loop plane cannot 
accurately be determined with the existing STEREO data, but they are roughly coplanar, 
based on the centroid motion constrained by AIA that is absent in EUVI/A from a near- 
perpendicular view (Fig. 5). Coplanar kink mode oscillations corresponds to a vertical 
polarization. 

2.5. Spatial Vatriation of Loop Oscillation 

In a next step we analyze the spatial variation of the transverse kink-mode oscillation 
a{t) as a function of the spatial loop position, which we specify with a segment number 
running from segment #1 at the loop length coordinate Si/L — 0.23 (near the first footpoint 
Fi) to segment #10 at S2/L — 0.78 (near the second loop footpoint). This analysis serves 
a two- fold purpose: (1) to detect possible asymmetries of the kink mode, and (2) to detect 
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possible propagating waves. 

In Fig. 6 we show the analysis of the loop oscillation of 10 different loop segments, 
numbered consecutively (# 1-10) from the southern loop footpoint Fi to the north-eastern 
footpoint F2 along the loop axis with loop length coordinate s. The location and orienta- 
tion of the time-slice stripes is indicated in the left panels of Fig. 6, the running-minimum 
difference time slices arc shown in the middle panels of Fig. 6, and the sinusoidal fits afit{t) 
to the loop amplitudes in the right-hand panels of Fig. 6, which also contains the best-fit 
parameters. If we discard the the two noisiest segments near footpoint Fi (Segment #1 and 
7^2 in Fig. 6), we obtain for the others a mean amplitude of (ai(s)) = 1.8 ±0.4 Mm, a mean 
period of {P{s)) = 373±30 s (6.2±0.5 min), and a mean starting time of {to{s)) = 399±35 s. 
Thus the variation of best-fit periods and starting times is only 8%, and thus we conclude 
that there is no significant phase shift of the oscillation amplitude along the loop that could 
be considered as a propagating wave. Thus, we deal with a pure standing wave of the fast 
MHD kink-mode. 

In Fig. 7 we show the spatial variation of the oscillation amplitude ai{s) in the context 
of an intensity image (Fig. 7, top left) and a running-minimum difference image (Fig. 7, 
top right). The locations of the 10 azimuthal time-slice stripes are shown in Fig. 7 (bottom 
left), over which the amplitude oscillation were measured in Fig. 6. The dependence of the 
oscillation amplitude ai(s) along the loop shows a maximum amplitude of a™"^ = 2.0 Mm 
near the loop apex. A sinusoidal displacement along the loop axis is ideally expected for 
a kink cigcn-modc with fixed nodes (compare with the analogy of a violin string). Our 
measurements, however, rather show a slightly distorted and asymmetric function, which 
can be approximated by a squared sine function (to account for the curvature of the loop) 
and a nonlinear dependence a(s°"^^) on the loop length (to account for the asymmetry of 
the loop, as evident from the stereoscopic triangulation of the footpoints, see Fig. 7 bottom 
right panel). 



where s = and s = L mark the nodes at the true footpoints Fi and F2. The observed oscil- 
lation amplitudes a(s) follow the squared sine function closely in the range of s/L < 0.6, but 
deviate in the range of 0.6 < s/L < 0.75, probably because of the interference of a secondary 
oscillating loop. 




(5) 
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2.6. Multiple-Loop Oscillations 

From watching the AIA 171 A movies (see electronic supplementary data to this paper) 
and from the time-slice plots shown in Figs. 2, 3, and 6 it appears that multiple loops 
are involved in kink-mode oscillations. The previous analysis has determined the average 
dynamic parameters of the collective ensemble of individual loop strands. In Figure 8 we 
are fitting a 2-loop model to the time-slice plots obtained near the loop apex, which yields 
slightly different periods (Pi = 375 s and P2 — 336 s), oscillation amplitudes (oi = 2.3 
Mm and 02 = 2.0 Mm), and loop centroid positions (xi = 6.8 Mm and X2 = 6.1 Mm), but 
a common start time = 423 s (i.e., 19:12:03 UT). Thus the loops are excited in phase, 
but the secondary loop has an oscillation period that is about 10% shorter. The secondary 
loop seems also to have a shorter lifetime and is only visible in the 171 A filter for about 2 
oscillation periods (compared with 3.6 periods of the primary loop.) 



2.7. Intensity Modulation During Loop Oscillations 

The 171 A intensity of the background-subtracted loop intensity exhibits strong modu- 
lations, being strongest near the beginning, but fading out gradually at the end of the time 
interval of oscillations. We show in Fig. 3 (bottom panel) the background-subtracted inten- 
sity flux profile fit) as measured near the loop apex from the Gaussian cross-sectional profile 
fits (Eq. 2). Amazingly, the intensity flux modulation appears to be in synchronization with 
the oscillation amplitude, which is a very interesting property that we have not noticed in 
previous observations of loop oscillations (e.g., in the 26 cases observed with TRACE; As- 
chwanden et al. 1999, 2002). In fact, the loop flux modulation f{t) (Fig. 3, bottom panel) 
occurs in anti-phase to the amplitude modulation a{t) that is measured in upward direction 
away from the loop curvature center. In addition to the oscillation-modulated variation, the 
flux decays as a function of time, which can be described with a linear decay rate df/dt, 
similar as found for 8 loops with kink- mode oscillations during the 2001 Apr 15, 21:58 UT, 
flare (Aschwanden and Terradas 2008). Thus, we flt a sinusoidal function with a linear decay 
rate /2/P/, 

'Mt-tf)\ , , (t-tf) 



ffuit) =fo + h sin [-^^^ ) + f^-rf^ ■ (6) 

We flnd a peak flux of /o = 61 DN s~^, a flux modulation of /i = 7.5 DN/s, a hnear decay rate 
of /2/P/ = -0.037 (DN s-2), which deflnes a loop lifetime of %e = ^/(/o//2) = 1649 s (27 
min) and is compatible to the loop cooling times Tcooi = 17±7 min found in Aschwanden and 
Terradas (2008), modeled also in Morton and Erdelyi (2009, 2010). It is therefore suggestive 
to interpret the observed lifetime of the oscillating loop as the detection time of a loop that 
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cools through the AIA 171 A passband. For the flux modulation that is anti-correlated with 
the amplitude oscillation we suggest an interpretation in terms of density compression by 
cross-sectional loop width oscillations, similar to a sausage mode, which is modeled in the 
next section. 



2.8. Density Modulation During Loop Oscillations 

In the previous Sections we established that the vertical oscillation amplitude amounts 
to a^"^/r;oop = 2/52 = 4%. If we assume that the loop is embedded in a magnetic field with 
a constant pressure across the loop cross-section (to first order) in a low plasma /9-parameter 
environment, the ambient magnetic field lines arc expected to oscillate in synchronization 
with a displacement that is proportional to the loop amplitude. A consequence of this 
scenario is that the loop cross-sectional radius r(t) — rO varies proportionally to the loop 
amplitude a{t), 

rit) =rJl + ^) , (7) 

leading to a modulation of the loop cross-sectional area A{t) — Trr^{t) that scales quadrati- 
cally to the loop radius r{t), 

A{t)^Ao(l + ^y . (8) 



^'loop , 

Since the loop footpoints are anchored at fixed positions in the photosphere, we can char- 
acterize the oscillating loop shape with an elliptical geometry to first order, which oscillates 
around the semi-circular geometry of the loop at rest, as depicted in Fig. 7. The loop length 
of a half ellipse is mathematically (to first order) , 

T f minor ~l~ major , /^\ 

L = 7r '— + ..., (9) 

where Tminor is the minor semi-axis and Vmajor is the major semi-axis of the ellipse. For the 
semi-circular limit the radii are equal, Tminor = t major = Tioop, and the half loop length is 
L = Tirioop- Assigning the minor axis to the half footpoint separation, r„,i„o, = ^ioop? and the 
major axis to the vertical radius with a small oscillation amplitude, r major = ''^'loop + ci{t), the 
eUiptical loop length varies (to first order) as. 



L{t) = TT {^loop + ^ + ...^ = Trr^oop {l + 0^ + 



(10) 



The volume of the loop, V{t), varies then consequently with the 5/2-power of the amphtude 
variation (to first order), 

V{t)=A{t)L{t) = vJl + ^ + y' . (11) 

\ floop J 
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The electron density inside the loop, assuming particle conservation in adiabatic compression 
and expansion processes, varies then reciprocally to the loop volume, 

ne{t)o,V{t)-^^no(l + ^ + ..) ' . (12) 

\ ^loop / 

For optically thin emission, as it is the case in EUV and soft X-rays for coronal conditions, 
the flux intensity scales with the square of the electron density times the column depth dz 
(which is here assumed to be proportional to the loop diameter dz oc -^yl(t)), yielding an 
anti-correlation of the flux with the 4-th power of the amplitude oscillation, 

F{t)o^nl{t)dz{t)^Fo(l + ^ + ..) . (13) 

\ ^loop / 

Thus, the small-amplitude variation of a^"'^ /rioop — 2.2/52 = 0.042 is amphfled with the 
4-th power, 

/ ^max \ 

Fma. = i^O 1 - 4^ + ... U 1.18-1 , (14) 

\ ^loop / 

which yields a flux modulation of 18% with respect to the mean value /q. In Fig. 3 (bottom 
panel) we fltted the flux variation and found indeed an average mean modulation factor of 
= /i/(/o(t)) ^ 7.45/30) ^ 24%, for the average of the total flux (/o(t)) ~ (61 + 10)/2 ^ 30 
DN s~i during the oscillatory episode. Thus our model predicts the correct time phase and 
approximate amount of oscillatory intensity flux modulation, which is anti-correlated to the 
sinusoidal loop amplitude oscillation (Fig. 9). The MHD wave mode that is associated with 
cross-sectional variation is called sausage mode or symmetric m = mode of fast MHD 
waves (e.g., Roberts 1984), which has a distinctly different eigen-mode period than the kink 
mode. The cross-sectional and density variation that were found in synchronization with the 
kink mode here (which has the same geometric and density properties as the sausage mode, 
but a different period than predicted by the MHD dispersion relation) , is a novel result of 
this study. This characteristic seems to be a particular property of oscillations in the loop 
plane (Fig. 7), also called "vertical polarization of kink mode" (Wang and Solanki 2004; 
Verwichte et al. 2006a,b), which would not occur (to flrst order) for transverse oscillations 
in perpendicular direction to the loop plane. 



2.9. Density and Temperature Analysis of Oscillating Loop 

Having the 6 coronal AIA fllters available that cover a temperature range of T f=:i 0.6 — 16 
MK for the entire oscillation episode with the same cadence of 12 s we are in an unprece- 
dented position to conduct an accurate diagnostics of the electron temperature and density 
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of the oscillating loop. For this purpose we extract loop-aligned subimages in all 6 coronal 
wavelengths in 10 spatial segments (sj,i = 1,...,10) at the loop locations as indicated in 
Fig. 7 (bottom left) and at 10 consecutive times {tj.j = 1, 10) during the time interval of 
19:05-19:35 UT. We show the 6 x 10 subimages for the 10 different times for loop segment ^ 
6 near the loop apex in Fig. 10 (left half), as well as the averaged cross-sectional loop profiles 
resulting from these subimages in Fig. 10 (right half). We calculate also the cross-correlation 
coefficients of these subimages with the simultaneous subimage in the detected wavelength 
of 171 A (indicated by the numbers in each subpanel in Fig. 10). From this information 
shown in Fig. 10 it is very clear that the oscillating loop exhibits a near-Gaussian cross- 
sectional profile only in the 171 A filter, while the 131 and 193 A filters show only a mild 
correlation [CCC ~ 0.4 — 0.5) and the remaining filters (211, 335, and 94 A) are absolutely 
uncorrelated {CCC ~ 0.0 — 0.1), which already narrows down the loop temperature to peak 
response temperature of the 171 A filter at T a; 10^'^ ~ 0.8 MK. 

In Fig. 11 (left side) we show the AIA temperature response functions, where the low- 
temperature response of the 94 A filter is corrected by an empirical factor of ^94 = 6.7 ± 1.7 
(Aschwanden and Boerner 2011). The total fluxes /*°*(t) (histograms with error bars in 
Fig. 11 middle panels) and background fluxes f^"-^^{t) (hatched areas in Fig. 11 middle panels) 
arc also shown, where the background is evaluated based on the Gaussian cross-sectional 
fits (Fig. 10). The difference f°"P{t) = f"\t) - /^"^^(t) is attributed to the EUV fiux of 
the oscillating loops and is modeled with a single-Gaussian differential emission measure 
(DEM) distribution by forward-fitting (according to the method described in Aschwanden 
and Boerner 2011), 



with the best-fit DEM solutions shown in Fig. 11 (top right panel) for the 10 consecutive time 
steps. The single- Gaussian DEM fits yield an average peak temperature of Tq = 0.57 ± 0.14 
MK and a Gaussian temperature width of <yiog{T) = 0.18 ± 0.10 (Fig. 11, right side), which 
corresponds to a near-isothermal temperature distribution at the limit of the temperature 
resolution Alog{T) 0.3 of the AIA filters, similarly as found for a statistical set of other 
loops analyzed from TRACE (Aschwanden and Nightingale 2005) or AIA (Aschwanden and 
Boerner 2011). The goodness-of-fit of the best-fit DEM solutions is found to be Xred — 
1.10 ± 0.05. The average agreement of the observed and modeled fiuxes is found to be < 3% 
in the three filters with the highest fiuxes (Fig. 11, middle column). The largest relative 
deviation occurs in the 94 A filter, which are known to have an incomplete temperature 
response function due to missing lines of Fc X transition (Aschwanden and Boerner 2011). 
Also we have to keep in mind that the largest flux deviations in the fits are only in the 
order of 0.2 — 0.5 DN/s in the 94, 131, and 335 A channels, which results mostly from 




(15) 
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uncertainties in the (time- variable) background evaluation rather than from the statistical 
photon noise. 

Assuming a filling factor of unity, we can estimate the mean electron density in the 
oscillating loop, 

for which we obtain a mean value of no = (1.9 ± 0.3) x 10^ cm~^, based on average loop 
widths oi w = 4.9 ± 0.6 Mm (Fig. 11, right side), which is measured near the apex of the 
loop for the segment # 6 shown as cross-sectional loop profiles in Fig. 10. 



2.10. Radiative Cooling Time Scale 

Since the issue has been raised that the lifetime of oscillating loops (defined by the 
detection time in a given temperature filter) is commensurable with the duration of an 
observed oscillation event (Aschwanden and Terradas 2008; Morton and Erdelyi 2009), let 
us explore whether the theoretically predicted time scales are consistent with the observed 
flux decay. At the relatively low coronal temperatures of < 1.0 MK observed in EUV, 
radiative cooling is the dominant time scale, while conductive cooling is only relevant at 
much hotter plasma temperatures in soft X-rays. Assuming an impulsive heating episode 
with subsequent cooling we can approximate the temperature evolution with an exponentially 
decaying function over some temperature range. 



Te{t - h) ^ Teih) exp 



{t-h) 



'^cool 

where the temperature cooling time TcooI corresponds to the radiative cooling time r^ad, 



(17) 



5 noAo 

1.0 ~ 10~^^'^^ erg cm^ s~^ being the radiative loss rate at EUV temperatures (T 1.0 
MK). For our measured values of To = 0.57 MK and no = 1.9 x 10^ cm~^ at the loop apex we 
estimate Trad ~ 2750 s (46 min). The loop lifetime Tin in the 171 A fllter, which has a FWHM 
temperature range of Tm = 0.53 — 1.17 MK, is then Tin = 2200 s (37 min) according to 
Eq. 17, which if fully consistent with the observed flux decay time Tufg = P/(/o//2) = 1650 
s (27 min) based on the fitted time profile (Eq. 6) to the measurement shown in Fig. 3 
(bottom panel). Therefore, we can interpret the observed fiux decay seen in the 171 A as a 
consequence of the radiative cooling time. Based on this cooling scenario we would predict 
an initial temperature of Te{t — ti) 0.8 MK at the beginning of the oscillation event and 
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a temperature of Te{t — t^) ~ 0.4 MK at the end of the oscillation episode. The predicted 
temperature drop Tf,{t) for a single loop cannot be retrieved by the DEM modeling (Fig. 11, 
left second panel), because the two nearly cospatial oscillating loops cannot be separated and 
thus we can only measure the combined emission measure-weighted temperature evolution of 
near-cospatial loops. Considering the analysis of the temperature evolution in other channels, 
only the 131 A channel has a temperature response to slightly lower values than the 171 A 
channel, but is about two orders of magnitude less sensitive and thus is unsuitable for a 
quantitative analysis. 



2.11. Excitation of Loop Oscillations 

The exciter or trigger of the loop oscillations is very likely the associated flare event to 
the east of the oscillating loop. If we calculate the projected distance between the flare site 
flare = iV flare — —410") (sectiou 2.2) and the apex of the oscillating loop {xapex = 

+698", yapex — — 243", /lapei; ~ 50") (sectiou 2.4), taking the stereoscopically triangulated 
3D loop position into account, we obtain and Euclidian distance of Lexc — 353" (256 Mm). 
Given the time delay between the start of the loop oscillation at the apex (19:12:12 UT) 
and the flare start (19:10:00 UT; ±6 s), we obtain the following 3D propagation speed of the 

Le.e 256,000 km ^q^..^.,, _i 

Vexc — tt:. — = = 1940 ± 125 km s , (19 

T,,^^ 132±6s ' ^ ' 

which is a typical Alfvenic (magneto-acoustic) speed in the solar corona. Thus, we can 
conclude that the loop oscillation is initially triggered by a fast MHD wave with Alfvenic 
speed. Moreover, the direction of the initial excitation in west-ward direction is in the same 
direction as the propagation direction of the fast MHD wave that propagates with Alfvenic 
speed concentrically away from the flare site. Although the angle between the Alfvenic wave 
direction and the loop oscillation amplitude cannot be determined with high accuracy, it is 
closer to parallel than perpendicular, as would be expected in a vortex-shedding scenario 
(Nakariakov et al. 2009), where the kink-mode oscillation occurs in perpendicular direction 
to the local plasma flow direction. In the case analyzed here, it appears that the plasma 
volume in the westward direction is stretched out in the same direction, probably following 
a narrow-angle cone of open magnetic held where the CME escapes. An associated CME 
on the south-west side of the Sun is visible in SOHO/LASCO and EIT movies. Generally, 
excitation of kink-mode oscillations with vertical polarization are rare, because they need 
special circumstances with an exciter near the curvature center of the loop (Selwa et al. 2011). 
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2.12. External/Internal Density Ratio of Oscillating Loop 

Coronal seismology allows us to determine the mean magnetic field in a loop in the 
kink-mode oscillation mode based on the internal Alfven speed va inside the oscillating 
loop, 

= -^L= . (20) 

which can be related to the (fundamental) kink-mode period Pkink using the phase speed Ck 
inside the flux tube (Roberts et al. 1984), 



_ 21/osc _ 2Losc H- + Pe/ Pi /r)-,\ 

Ck V 2 

which depends on the total length Lose of the oscillating loop and the densities externally 
(ne) and internally (n^) of the loop. 

On the other hand, we can estimate the external Alfven speed Vjye, which depends on 
the external magnetic field Be and density Ue, 

VAe = . (22) 

If we interpret the exciter speed v^xc between the flare site and the (apex) location of the 
oscillating loop (Eq. 19) as an Alfvenic wave, we obtain a direct measurement of the external 
Alfven speed (supposed the wave is not super- Alfvenic) , 

Lexc 

VAe = Vexc = Tf; ■ (2<3j 

exc 

Moreover, since the magnetic pressure is generally dominant over thermal pressure in the 
solar corona, the magnetic field Bi internally and B^ externally to the loop boundary have to 
match for a self-consistent magnetic field model of a loop embedded into an external plasma. 
Based on the definitions of the Alfven speeds (Eqs. 20, 22), the ratio of the Alfven speeds 
depends then only on the density ratio, 

Be^B,^^ = M, (24) 

Va V 

which can be directly determined from the kink-mode period, loop length, and exciter speed 
with Eqs. (21) and (23), 



'^i 1 / L^xc Pkink 



Tip 2 \ Lqsc ^exc 



1 . (25) 



- 15 - 



Some uncertainty arises from the unknown positions of the actual nodes of the kink mode. 
Since the oscillation has only been observed over a loop segment of Lg^g = 123 Mm, which 
is a lower limit of the node separation, while the stereoscopically triangulated full loop 
length Lioop down to the solar surface represent an upper limit, we might estimate a realistic 
value with uncertainty from the arithmetic mean of the two limits, i.e.. Lose ~ {Lioop + 
Lseg)/2 = 143 ± 20 Mm. Thus, for the physical parameters determined above {Lexc = 256 
Mm, Tea;c = 132 ± 10 s), and the best-fit values of the primary oscillating loop (Fig. 8: 
Lose ~ {Lioop + Lgeg)/'^ — 143 ± 20 Mm, Pkink — 375.6 s), we obtain a value of rij/rie ~ 
12.0 ± 1.8, or an inverse ratio of ne/rii 0.08 ± 0.01. This value is commensurable with 
alternative methods, using stratified hydrostatic density models of the background corona, 
for which a statistical average of n^jni = 0.30 ± 0.16 was found. (Aschwanden et al. 2003). 
Consequently, for a density ratio of rie/rii ^ 1/12, the ratio of the Alfven speeds is then 
expected to be (vAe/vA) = ~ ~ 3.4. In our case, the external Alfven speed is 

vac ~ 1940 ± 100 km s~^, and the internal Alfven speed is va ~ 560 ± 100 km s~^. 



With this novel method of measuring the density ratio rie/ni from Alfvenic propagation 
speeds external and internal to the fiuxtube, the magnetic field in the oscillating flux tube 
and immediate surroundings is then fully constrained with Eqs. 20 and 21 (Nakariakov et 
al. 1999), 



for which we obtain S = 4.0 ± 0.7 G, based on the measurements of Lose — (143 ± 20) x 10^ 
cm, Pkink = 375.6 s, rij = (1.9 ±0.3) x 10^ cm~^, and the density ratio (ne/rii) ^ 0.08 ±0.01. 

This magnetic field scenario can be tested with observed magnetic field data from the 
Helioseismic and Magnetic Imager (HMI) on SDO. In Fig. 12 (bottom) we show a HMI 
magnetogram recorded at 19:04:16 UT, at the beginning of the analyzed time interval. The 
fiare location is situated in the core of the AR, right at the neutral line with the largest 
magnetic fiux gradient, while the oscillating loop is located beyond the western boundary 
of the active region in a low magnetic-field region that is governed by a "salt-and-pepper 
pattern" of positive and negative magnetic pores (see enlargement in Fig. 13 top left). A 
potential-field source surface (PFSS) model calculation is shown in Fig. 12 (top panel), which 
is dominated by a bipolar arcade above the neutral line in east-west direction. 

The magnetic field in the environment of the oscillating loop can be modeled with 
potential-field or non-potential field models, but both are known to show misalignments 



2.13. Magnetic Field Modeling 




kink 



(26) 
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with the 3D geometry of stereoscopically triangulated loops of the order of amis ~ 20° — 40° 
(DeRosa et al. 2009; Sandman et al. 2009), while simple potential field models calcu- 
lated from a small set of unipolar magnetic charges (Aschwanden and Sandman 2010) or 
magnetic dipoles (Sandman and Aschwanden 2010) achieved a reduced misalignment of 
amis = 13° — 20°. For a simple plausibility test of the magnetic field strength inferred from 
coronal seismology, we model the 3D field at the location of the oscillating loop with an an- 
alytical model of two unipolar charges with opposite magnetic polarities that are buried in 
depths Zi and Z2 and have maximum longitudinal magnetic field strengths of Syi = -1-187 G 
and Bp — —63 G at the observed positions {xi, yi) and {x2, ^2) of the nearest magnetic pores 
in the HMI magnetogram (marked with circles in Fig. 13, top left panel), where {x, y, z) is a 
cartesian coordinate system with the xy-plane parallel to the solar surface. The correspond- 
ing absolute field strengths vertically above the buried charges are Bi = -Byi/ cos('i?i) = 296 
G and B2 = Bp/ cos('i?2) = —89 G, where 'dj are the line-of-sight angles. Thus, in this model 
we have only the two free variables of the depths Zi and Z2 to fit the model of resulting mag- 
netic field lines to the observed loop. The magnetic field resulting from the superposition of 
two unipolar magnetic charges is then given by (Aschwanden and Sandman 2010), 



in terms of the vector = [{x — xj), {y — yj), {z — Zj)], with Xj = {xj,yj,Zj) being the 
locations of the buried unipolar magnetic charges and Bj the magnetic field strength at 
the solar surface above the magnetic charges. The ratio of the two free variables Zi and 
Z2 determine the asymmetry of the field lines. For the observed oscillating loop we find 
values of Zi = 0.5" and Z2 = 1.5" to reproduce approximately the observed shape (Fig. 13, 
bottom left). The field line that closest fits the projected location of the oscillating loops 
has magnetic field strcnghts of Bi = 296 G and i?2 = 29 G at the photospheric field line 
footpoints and i? = 6 G at the apex, which compares favorably with the magnetic field 
strength of B^ink = 4.0 ± 0.7 G deduced from coronal seismology. 

However, since the magnetic field B{s) varies along the loop, the Alfven speed varies 
proportionally and the Alfvenic transit time during one oscillation period is given by 



which defines an average magnetic field {B) that is equivalent to a fluxtube with the same 
kink-mode period P and a constant magnetic field value {B) by 




(27) 




(28) 




(29) 
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for which we obtain (B) — 11 G, which is a factor of 1.8 higher than the minimum value at the 

apex, Bapex = 6 G, or a factor of 2.8 higher than inferred from seismology, -Bfcmfc = 4.0 ± 0.7 
G. This difference between the seismological and magnetogram-constrained magnetic field 
value, derived for the first time for an oscillating loop here to our knowledge, is perhaps not 
too surprising, given the ambiguity of potential-field models, non-potential field models, and 
uncertainties of the footpoint locations (which require stereoscopic information). 

3. DISCUSSION 

In this well-observed loop oscillation event, which we analyzed with AIA/SDO, HMI/SDO, 
and EUVI/ STEREO, we derived a comprehensive number of physical parameters (listed in 
Table 1) that could not be determined to such a degree in previous observations. In the fol- 
lowing discussion we compare the observational results with theoretical models, predictions, 
and discuss interpretational issues. 

3.1. Coupled Kink-Mode and Cross-Sectional Oscillations 

The basic theory for fast magneto-acoustic waves, which predicts kink and sausage eigen- 
modes for slow (acoustic) and fast (Alfenic) MHD waves, has been derived for a straight 
(slender) cylindrical fiuxtube (e.g., Edwin and Roberts 1983). For such an idealized geom- 
etry, the periods of the fast kink and sausage mode have quite different regimes, and the 
sausage mode has a wavenumber cutoff with no solution of the dispersion relation for ka <1 
(with k the wave number and a the fiuxtube radius), which corresponds to a cutoff at a 
phase speed of Vph = vab- From this theory, no sausage eigen-mode is predicted for periods 
that correspond to kink-mode oscillations, Pfcmfc = 2L/v^. In contrast, our analysis clearly 
demonstrates the presence of a kink-mode with coupled sausage-like behavior, as measured 
by the cross-sectional loop width variations and anti-correlated density variations. The ques- 
tion arises why this dynamical behavior is not predicted by existing theory? One possible 
explanation is that the loop length is not constant but changes as a function of time in syn- 
chronization with the transverse oscillation amplitude. This is most plausibly seen in Fig. 1, 
where the excitation direction originating from the flare location propagates in approxi- 
mately the same direction as the loop plane, and thus excites a significant component of the 
"vertical" polarization mode (i.e., the loop plane and the oscillation plane are near-parallel), 
as inferred for one case in Wang and Solanki (2004) and analytically studied in Verwichte 
et al. (2006a,b). Most kink-mode oscillations have a horizontal polarization, as determined 
with STEREO (e.g., Verwichte et al. 2009), but density oscillations have alse been noted 
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in previous kink-mode oscillations (e.g., Verwichte et al. 2009, 2010). If the loop oscillates 
in vertical polarization, the length of the loop can vary during the kink-mode oscillations, 
with a linear dependence on the oscillation amplitude to first order (in the elliptical approx- 
imation, see Eqs. 9 and 10). Thus, the periodic shrinking and stretching of the fiuxtube is 
likely to cause a bulging and thinning of the central loop cross-section, which is exactly what 
a sausage mode does. A consequence of the length variation L{t) is also a magnetic field 
variation B{t), which scales reciprocally to the cross-sectional area of the sausage mode, i.e., 
B{t) oc A~^(t), due to the conservation of the magnetic flux, i.e., B(t)A(t)— constant. 

The coupling of kink-mode and (sausage-like) cross-section and density oscillations thus 
might be a special case that occurs only when the loop length is varied, which most likely 
occurs for vertical polarization and requires an initial excitation in direction of the loop 
plane. It would be interesting to investigate this prediction of coupled cross-section and 
density oscillations as a function of the exciter direction or kink-mode polarization, which 
depends on the location and orientation of the loop plane with respect to the propagation 
direction of a flare or CME-related disturbance. Since CME bubbles and erupting flux ropes 
get stretched out during the initial expansion, it is natural that ambient magnetic field 
lines become stretched too, which applies also to oscillating loops. Statistics on different 
polarization types of kink-mode oscillations is still small and their identification based on 
difference images is often ambiguous (Wang et al. 2008). 

An alternative interpretation of the amplitude-correlated flux variation is an aspect- 
angle change of the oscillating loop, which causes a variable line-of-sight column depth of 
the loop diameter w(t) = wocos['d{t)] (Cooper et al. 2003), and hence would introduce 
a variation of the optically-thin EUV flux f{t) oc nl{t)w{t). However, the observed flux 
variation with a mean of 24% would require an aspect angle change of Ai? 40°, which 
is inconsistent with the observed stationarity of the loop shape during the entire oscillation 
episode. 

3.2. Multi-Loop Oscillations 

Evidence that multiple loops or strands are involved in this oscillation event is shown in 
Fig. 8, where we found shghtly different periods (by ^ 10%), amplitudes, centroid positions, 
and possibly different lengths (although not directly measured). The eigen-modes in a two- 
slab system was studied in ArrcgTii ct al. (2008) and it was found that the kink-mode periods 
may differ from a single loop when the distance between the loops is less than a few loop 
diameters. In our case, the projected centroid position is displaced by Ax = 0.7 Mm, 
while the loop diameters are w ~ 4.9 ± 0.6 Mm, so they could be close to each other. 
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Luna et al. (2008) simulated numerically the MHD behavior of two parallel loops and found 
four collective modes, kink (asymmetric) and sausage (symmetric) modes in both parallel 
and perpendicular direction to the plane that contains the axis of both loops, with four 
different frequencies, which is a generalization of the two modes of a single-loop oscillation. 
However, analytical solutions of a two-loop system yields only two different frequencies (Van 
Doorsselaere et al. 2008), which might differ from the numerical results of four different 
frequencies (Luna et al. 2008) due to the neglect of higher-order terms (Ruderman and 
Erdelyi 2009). 

A multi-threaded model with four loop threads was modeled with a 3D MHD code 
(Ofman 2009). For parallel threads, the evolution of the ensemble exhibits the same period 
and damping rate as a single loop, but for twisted threads, the periods become irregular and 
the damping much stronger, which seems not to apply to our case here. Either the multiple 
loops are near-parallel or sufficiently distant to each other. 

Resonant absorption in comphcated multi-strand loops was investigated by Terradas et 
al. (2008) and it was found that the damping behavior is not compromised by the complicated 
geometry of composite loops. One theoretical prediction of multi-loop oscillations is that the 
collective width w{t) increases with time due to a shear instability (Terradas 2009), but we 
do not observe such an effect (Fig. 11, bottom right panel), either because the two oscillating 
loops are not in sufficiently close spatial proximity or because the lifetime of the oscillating 
loops in the detected wavelength is too short. 



3.3. Damping by Resonant Absorption 

An unusual property of this oscillation event is that we do not observe any significant 
damping of the kink-mode amplitude over the duration of the oscillatory episode, so the 
ratio of the damping time to the period must be much longer than the observed number of 
periods, i.e., t^/P ^ A. This is in contrast to a statistical sample of 11 well-observed events 
with TRACE, where strong damping was found to be the rule, i.e., with to/P ~ 1-8 ± 0.8 
(Aschwanden et al. 2002). 

Resonant absorption as a damping mechanism for kink-mode oscillations was considered 
in Goossens et al. (2002). The ratio of the damping time to the period was calculated for 
resonant absorption by Ruderman and Roberts (2002) for a thin-boundary layer and by Van 
Doorsselaere et al. (2004) for thick boundaries. 
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where qtb ~ 0.75 is the correction factor for the thick-boundary layer, Igkin is the skin depth 
or thickness of the loop boundary that contains a density gradient, and g„ = ng/rij is the ratio 
of the external to the internal electron density in the loop. This density ratio was previously 
measured to g„ = 0.30 ± 0.16, based on loop flux intensities and hydrostatic models of the 
background corona (Aschwanden et al. 2003), and a skin depth ratio of rioop/lskin = 1.5 ±0.2 
was inferred, and hence the typical ratio of the damping time to the oscillation period was 
found to be r^/P fa 1.3. 

In our case, a similar density ratio of g„ ^ 0.08 was measured. Wc can reconcile 
the observed long damping time ratio of r/j/P ^ 4 only with a very small skin depth of 
hkin/f^ioop ^1/4. While previously analyzed kink- mode oscillations with TRACE exhibited 
typical temperatures of Tg ~ 1.0 — 1.5 MK, we deal here with a significantly cooler loop with 
a temperature of Tg fa 0.5 MK. It appears that such cooler loops have either a smaller skin 
depth or larger loop diameters than the warmer coronal loops, but no hydrodynamic model 
is known that predicts such an effect. 

3.4. Magnetic Field Comparisons 

Coronal seismology determines the magnetic field strength by setting the kink-mode 
period Pkink equal to the Alfvenic crossing time 2L/va forth and back along the loop length 
L, which yields a relationship for the magnetic field Bkink as a function of the loop length L, 
period Pkink, internal Ui and external density Ue (Eqs. 21 and 26). This method is one of the 
foundations of coronal seismology, initially applied by Roberts et al. (1984), Aschwanden et 
al. (1999), and Nakariakov and Ofman (2001). In principle, this analytical relationship can 
be put to the test by 3D MHD simulations of kink-mode oscillations of a plasma fluxtube by 
comparing the theoretical with the experimental values of the kink-mode oscillation periods 
Pkink or magnetic fields B. Such a test was conducted by DeMoortel and Pascoe (2009), but 
surprisingly the coronal seismology formula predicted a field strength {Bkink = 15 — 30 G) 
that was about a factor of 1.5 higher than the input values of P = 10 — 20 G of the MHD 
simulation. 

Here we attempted to vahdate the seismological magnetic field value {Bkink = 4.0 ± 0.7 
G) by a potential-field model that consists of two unipolar magnetic charges with opposite 
polarities that are buried near the footpoints of the oscillating loop and are constrained by 
the longitudinal magnetic field strengths observed in HMI magnetograms. The best-fit field 
line yielded a magnetic field value of Bapex = 6 G, which is a factor 1.4 higher than the 
seismological value. If we correct for the variable Alfven speed along the loop, we predict 
a seismological value of Bavg = 11 G, which is a factor of 2.8 higher than the theoretical 
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value. We note that the discrepancy of the our best-fit potential field model is in opposite 
direction to the discrepancy found from 3D MHD simulations by DeMoortel and Pascoe 
(2009). We believe that the discrepancy from magnetic field modeling methods mostly stems 
from the uncertainty of the footpoint locations, the spatial resolution of magnetograms, and 
the ambiguity of potential and non-potential field models. Stereoscopically triangulated 
loop oscillations hold the promise for obtaining more accurate measurements of the loop 
length and footpoint location. The most powerful self-consistency test needs to employ a 
combination of stereoscopy, numerical 3D MHD simulations, coronal seismology theory, and 
analytical magnetic field models. 



4. CONCLUSIONS 

Here we present the first analysis of a loop oscillation event observed with AlA/SDO, 
which occurred on 2010-Oct-16, 19:05-19:30 UT. The capabilities of AIA enable us for the 
first time to study such an event with sufficiently high cadence, spatial resolution, and 
comprehensive temperature coverage, which enables us to derive all important physical pa- 
rameters. In addition, magnetic modeling with HMI data can validate the magnetic field 
measurements based on coronal seismology. The major observational findings, interpreta- 
tions, and conclusions are: 

1. A flare with an associated CME that escapes the Sun along a narrow cone (in westward 
direction) excites kink-mode oscillations with a period of P = 6.3 min in a loop at a 
distance of Lexc = 256 Mm away from the flare site, after a time delay of Texc = 132 s, 
which yields an exciter speed of Vexc = Lexc/Texc ~ 1900 km s^^, which we interpret as 
a magneto-acoustic wave with Alfvenic speed and can be used as a direct measurement 
of the average external Alfven speed vac — v^xc outside the oscillating loop. 

2. The direction of the excitation and kink-mode oscillation amplitude is about in the 
same direction as the loop plane, which corresponds to a vertical polarization of the 
kink mode, causing a periodic stretching of the loop length and coupled cross-section 
and density oscillations, evident from the compression and rarefaction of the density, 
which produces an intensity variation that is amplified with the fourth power of the 
amplitude displacement. This behavior of kink modes with coupled cross-sectional and 
density variations are unusual and perhaps occur only in vertical polarization. They 
are not predicted by theory, which needs to be generalized for temporal variations of 
the loop length L{t). 
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3. There is evidence for a multi-loop system that is involved in the coupled kink and 

cross-sectional oscillations, consisting of at least two loop strands that have slightly 
different periods (~ 10%) but are excited in phase at the beginning. The fact that the 
two major oscillating loop strands arc not synchronized to the same period indicates 
a spatial separation of more than a few loop diameters. 

4. A full DEM analysis with all 6 coronal AIA temperature filters yields a temperature 
of T 0.5 MK and a density of ng ~ 2 x 10^ cm~^. Consequently, the loop oscillations 
are primarily observable in the 171 A filter, very faint in the 131 and 193 A filter, 
and essentially undetectable in the other filters. From this temperature and density 
measurement we estimate a radiative cooling time of Trad = 46 min, which explains the 
loop lifetime of Tuf^ = 27 min in the 171 A filter. 

5. The measurement of the external Alfven speed vac ~ 1900 km s~^ from the exciter 
speed and the internal Alfven speed = 560 km s^^ from the kink-mode period 
provides a direct measurement of the density ratio external and internal to the loop, 
ne/rii = 0.08 ± 0.01, which is commensurable with earlier hydrostatic models of the 
background corona {ue/ni — 0.30 ± 0.16; Aschwanden et al. 2003). This value pro- 
vides a fully constrained magnetic field measurement of the oscillating loop by coronal 
seismology, B^nk = 4.0 ± 0.7 G. 

6. For an independent estimate of the magnetic field in the oscillating loop we used a 
potential-field model with two unipolar magnetic charges, constrained by the photo- 
spheric magnetic field strengths {B^ = +296 G, B2 = -89 G) obtained from HMI/SDO 
magnetograms near the footpoints of the oscillating loop, which were localized by 
stereoscopic triangulation from STEREO/EUVI-A images. A best-fit model yields a 
magnetic field strength of Bapex = 6 G at the loop apex, or Bavg = 11 G when averaged 
along the loop. This independent test validates the coronal seismological value within 
a factor of ai 2. 

7. The oscillating loop exhibits no detectable damping over the observed four periods, 
which is unusual, compared with the statistical values of td/P = 1.8 ± 0.8 found from 
previous measurements. Damping by resonant absorption can only be reconciled with 
this observation if the skin layer (of the density gradient at the loop boundary) is much 
smaller than the loop radius. It is not clear if this property is a consequence of the 
unusual low loop temperature oiT ^ 0.5 MK. 

The excellent quality of the AIA data have provided more physical parameters of a 
coronal loop oscillation event than it was possible to determine in previous TRACE obser- 
vations, especially due to the much better cadence of 12 s, which allows us also to resolve 
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multi-loop oscillations, spatially and temporally. The measurements of more physical pa- 
rameters provide stronger constraints on the theory and raise new problems that need to 
be addressed by analytical theory or MHD simulations: (1) What is the 3D geometry and 
timing of the exciter mechanism and how does it affect the polarization of kink-modc os- 
cillations? (2) Can we explain the coupling of kink mode and (sausage-like) cross-sectional 
and density oscillations? (3) Can we explain kink-mode oscillations with no damping? (4) 
How do multi-loop oscillations interact with each other and how do the MHD wave modes 
couple? (5) How accurate are magnetic field measurements based on coronal seismology and 
how can they be validated with magnetic field models? Progress in these questions calls 
for modehng that combines stereoscopy, numerical 3D MHD simulations, coronal seismology 
theory, and analytical magnetic field models. 
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Table 1. Observables and physical parameters of analyzed loop oscillation event 



Paraiiieter 



Date of observations 

Time interval of analyzed observations 

Time range of GOES flare 

Flare onset of impulsive phase 

Start of loop oscillations 

GOES flare class 

Active region number 

Flare location 

Location of oscillating loop footpoints 

Location of loop apex 

Distance of flare to loop apex L^xc 

Delay of flare start to loop oscillation Tg^c 

Exciter speed Vexc — VAe 

Height of loop apex hapex 

Distance from Sun center 

Full loop length Lioop 

Length of oscillating loop segment Lseg 

Loop curvature radius r/oop 

Loop FWHM diameter w 

Loop inclination angle to vertical 

Polarization angle of kink oscillation 

Drift velocity of loop centroid ds/dt 

Oscillation period of loop P 

Oscillation amplitude of loop Oi 

Number of oscillation periods Np 

Loop lifetime r^/e = fo/idf/dt) 

Ratio of loop amplitude to radius d'l'"'''^ /rioop 

Observed flux modulation /1//0) 

Electron temperature Tg 

Temperature width (Tiog{T) 

Electron density rig 



2010-Oct-16 
19:05-19:35 UT 
19:07-19:12 UT 
19:10:00 (±6 s) UT 
19:12:12 (±6 s) UT 
M2.9 

NOAA 1112 
[390",-410"], W26/S20 
[685",-305"], [615",-268"] 
[698",-243"] 
275 Mm 
132 ± 10 s 
1940 ± 100 km s^^ 
37 Mm 
740" (O.77i?0) 
163 Mm 
123 Mm 
52 Mm 
4.9 ± 0.6 Mm 
20° ± 20° 
vertical 
0.8 km/s (towards west) 
375.6 s (6.3 min) 
1.7 ±0.4 Mm 
3.6 

1650 s (27 min) 
0.042 

0.24 (0.18 predicted) 
0.57 ±0.14 MK 
0.18 ±0.10 

(1.9 ±0.3) X 10^ cm-3 
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Table 1 — Continued 



Parameter 


Value 


External Alfven speed Vexc = vac 


1940 ± 100 km s"^ 


Internal Alfven speed va 


560 ± 100 km s"^ 


External/internal density ratio Ue/rii 


0.08 ±0.01 


Magnetic field at loop apex B^pex 


4.0 ±0.7 G 


Magnetic field at loop footpoints Bfoot 


+296, -89 G 


Damping time ratio Tdamp/ P 


>4 



200 400 600 800 

EW (arcsec from Sun center) 




EW (arcsec from Sun center) 

Fig. 1.— AIA 171 A image of flare observed on 2010-Oct-16 19:22:36 UT shown with the flux 
on a logarithmic flux scale (top panel) and as difference image with respect to 19:21:00 UT 
(bottom panel). The flare location is marked with a cross (in the center of the diffraction 
pattern) and a box indicates the location of the oscillating loop. [See also movies in 171 
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File = 201 0-Oct-1 6_1 905_1 71_a 

a) subtract boxcar-smoothed image (nsm= 5) 




500 1000 1500 

c) one-sided running difference, F(ti)-F(ti.i), j= 5 




500 1000 1500 

Time t(s) 



Fig. 2. — Five different enhancement algorithms to visuahze the loop oscillation in a time- 
slice plot with the time running along the x-axis (time range is 2010-Oct-16 19:05 — 19:35 
UT) and the spatial coordinate (defined in Fig. 1 along the stripe perpendicular to the loop 
through the loop apex): (a) highpass filtering by subtraction of a boxcar-smoothed image 
(top panel); (b) baseline difference by subtraction of the first time slice (second panel); 
(c) one-sided running time difference (third panel); (d) symmetric running time difference 
(fourth panel); and (e) minimum running time difference (bottom panel; At = ±5 frames). 
Thftrp annfiar snmp 'Whnfis" nr "miiltinlft nfirinds" ft pr arnimd t Ri 1000 
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Fig. 3. — Time-slice diagram of 171 A flux of oscillating loop with the running-minimum 
difference method (top panel; At = ±10 frames), with Gaussian cross-sectional fits (second 
panel) and best-fit model (third panel), based on a fit of the transverse oscillation amplitude 
a{t) with a sine function plus linear motion (fourth panel), and anti-correlated flux modu- 
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Distance x [Mm] 



Fig. 4. — Cross-sectional flux profiles F{x,t) (blue) obtained from the running-minimum 
difference technique (Fig. 3, top panel) and Gaussian fits (red). 
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Fig. 5. — SDO and STEREO observations of oscillating loop: An AIA 171 A differ- 
ence image (19:19:24-19:21:00 UT) is shown at the bottom left, and a near-simultaneous 
STEREO/EUVI/A 195 A difference image (19:15:30-19:20:30 UT) with additional highpass 
filtering is shown at the bottom right. The corresponding field-of- views and loop outlines 
are shown for both instruments in the top panels. The thick solid curve in the AIA image 
indicates the tracing of the oscillating loop segment, which is fitted to the corresponding 
Innn seerment in ETIVT/A bv 3D cnnrdinate tran.sfnrmatinns with variable altitudes and in- 
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Fig. 6. — Ten different stripes perpendicular to the loop spine (left panels), used to extract a 
running- minimum difference plot (middle panels), and sinusoidal fits of the loop oscillation 
amplitude as a function of the loop position (right panels), from loop segment 1 (near 
footpoint Fi) to loop segment 10 (near footpoint Fo). 
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Fig. 7. — Location of stereoscopically triangulated loop footpoints (marked with diamonds 
and labeled with Fi and F2) and 10 loop segments with the rectangular subimages over which 
the 10 time slices (shown in Fig. 6) were extracted (bottom left panel). The corresponding 
AIA 171 A image (top left) and difference image 19:21:00-19:19:24 UT (top right) are also 
shown. The magnitude of the transverse kink-mode oscillation amplitude is indicated with 
thick bars (bottom panels), which fit a stretched sine function (Eq. 5). 
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Fig. 8. — Two-loop oscillation model fitted to the same oscillation amplitudes a{t) as mea- 
sured in Segment #6 shown in Figs. 3 and 4. The amplitudes of the two oscillating loops are 
indicated with solid and dashed curves (bottom panel). The two- loop solution is also visu- 
alized as time-slice plots for the absolute flux (third panel) and running-minimum difference 
representation (second panel). 
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Fig. 9. — Schematic of sinusoidal modulation of the loop vertical radius r{t) = rioop + ci{t), 
volume V{t) — Vq, density ne(t) — no, and flux intensity I(t) — Iq during one full oscillation 
period. Note that the density and intensity vary in anti-correlation to the loop amplitude, 
indicated with the grey-scale shading of the loops. 
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Fig. 10. — Cross-sectional loop profiles in 6 AIA filters and at 10 different times during the 
oscillation episode from 19:05 to 19:35 UT. The wavelength of the primary loop detection is 
171 A, with which the subimages are cross-scorrelated (with the cross-correlation coefficients 
given in each panel). 
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Fig. 11. — DEM modeling at 10 different times during the oscillation episode from 19:05 
to 19:35 UT: AIA response functions (left), flux versus time (middle: histogram), with loop 
background (middle hatched) and best-fit fiuxes (diamonds), DEM for 10 times (right top), 
temperature T(t) (right second panel), temperature width arit) (right third panel), electron 
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Fig. 12. — HMI magnetogram of same ficld-of-view as shown in Fig. 1, with potential field 
source surface (PFSS) model field lines (top) and locations of flare (diagonal cross), oscillating 
loop segment (curve), and stereoscopically triangulated footpoints (crosses) indicated. 
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Fig. 13. — Enlarged field-of-view (identical to Figs. 5 and 7) of the HMI magnetogram (top 
left), AIA 171 A image (top right), and HMI-based dipolar potential field model (bottom 
right) of oscillating loop (white curves). A field line that closely coincides with the oscillat- 
ing loop is shown separately (bottom right; black curves), constrained by the longitudinal 
magnetic field observed in the HMI magnetograms with — 187 G and Bp — —63 G. 



